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Abstract 


The mixing performance in a mixing section of a rich 
burn/quick mix/lean burn (RQL) combustor has been calcu- 
lated using a 3-D numerical model in a non-reacting environ- 
ment. The numerically calculated results have been com- 
pared with the measured data reported by Hatch, Sowa, 
Samuelsen, and Holdeman, 1992. The numerical 3-D tem- 
perature fields qualitatively agree with the experimental data. 
Also the development of the mixing flow and temperature 
non-uniformity trends throughout the mixing section for the 
numerically calculated results quantitatively agree with the 
measured data. The numerical model predicts less mixing 
and enhances the temperature gradients as compared to the 
measured data for the cases reported by Hatch et al. (1992) 
which include circular and slot orifice shapes (with different 
slant angles and aspect ratios). The predicted and measured 
results generally agree in the selection of the slanted slot 
orifice configuration yielding the best overall mixing perfor- 
mance (based on temperature uniformity) of all the configu- 
rations analyzed herein. 

Nomenclature 


A m = duct crossectional area, also A tot , m 2 
ACd = effective orifice area, m 2 

AR = area ratio (jet/mainstream) = Aj/A m = ACd/A m 
AMIX = area weighted overall temperature deviation 
from Teq, Eq 3 

AHOT = area weighted temperature deviation above 
Teq, Eq 4 

ACOLD = area weighted temperature deviation below 
Teq, Eq 5 

DR = density ratio (jet/mainstream) 
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DP/P 

f 

j 

J 

k 
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= total pressure loss across the mixing wall 
= non-dimensional temperature, Eq 2 
= radial vector direction 
= momentum-flux ratio (jet/mainstream) = 
M 2 /DR. also J =(MR)2/[(DR)( A/A m )2| 

= tangential vector direction 
= slot length (long dimension), m 


m tot = overall mixer mass flowrate, kg/sec 

= +m ma„, 
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mass flux ratio (jet/mainstream) = DR V:, t 
/U J 

' '-'main 

mass flow weighted overall temperature devia 
tion from Teq, Eq 6 

mass flow weighted temperature deviation 
above Teq, Eq-7 

mass flow weighted temperature deviation 
below Teq, Eq-8 

mass flowrate ratio (jet/mainstream) 
total pressure, atm. 

radial interface between the counter-swirling 

flows. Figure- 11 

radius of the mixing section, m 

radial distance from the centerline of the mixer, 

m 

temperature, K 

local axial velocity, m/sec 

approach mainstream axial velocity, in/sec 

velocity ratio (jet/mainstream) =Vjet/U main 

radial velocity of the jet, m/sec 

slot width (short dimension), m 

axial distance from the leading edge of the 

orifice , m 


p = fluid density, kg/m 3 

<(> = slot slant angle with respect to the axial 

direction, degrees 

Subscripts 
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Introduction 

Advanced highly fuel-efficient commercial turbo propulsion 
gas turbines pose a number of design challenges for combus- 
tion system including durability, aerothermal performance, 
wide operability range, and exhaust emissions. Gas turbine 
combustion systems in current use employ a single stage 
combustion zone which has been optimized to produce low 
exhaust emissions of carbon monoxide, and unburned hydro- 
carbons at idle and near-idle operating conditions. Exhaust 
smoke below the threshold of plume visibility throughout the 
engine operation range has also been maintained without 
adversely affecting engine starting, stability, and relight 
requirement. Moderate reductions in NO x emissions (up to 
30%) have been achieved in single-combustion- zone com- 
bustors by reducing unmixcdness and combustion zone 
residence time. Further reductions in high-power NO x 
emissions (e.g. up to 50% from state-of-the-art levels) would 
require design and development of two-stage combustion 
concepts, namely, radially or axially staged combustion. 
These types of combustion systems are currently under 
development. 

In order to achieve more demanding (70 to 90%) NO x reduc- 
tion goals of the advanced turbine engines which will have 
considerably higher overall pressure ratios and turbine rotor 
inlet temperatures, significant advances are needed in the 
design and development of combustion systems employing 
multi-staging concepts including lean premix (LPP) and rich 
burn/quick mix/lean burn (RQL) (Shaw, 1991). The RQL 
concept studied in-depth by Rizk and Mongia (e.g. Rizk and 
Mongia. 1990; 1991; 1993) appears to have a number of 
advantages over other competing ultra-low NO x concepts 
provided we can significantly reduce the NO x formation by 
enhancing the mixing in the transition from rich to lean-side 
combustion. An experimental and analytical research effort 
under the sponsorship of NASA Lewis Research Center (e.g. 
Smith. Talpallikar. and Holdeman. 1991: Talpallikar et al., 
1991; Vranos et al.. 1991; Hatch et al., 1992, Bain, Smith, 
and Holdeman, 1992; Liscinsky et at.. 1992; ; Oechsle, 
Mongia, and Holdeman. 1992; Kroll et al., 1993; Liscinsky. 
Vranos, and Lohmann, 1993) is underway to study and 
identify the critical design and flow parameters affecting the 
mixing effectiveness. 

The 3-D analytical technique can be an advantageous tool 
since a large number of configurations can be easily modeled 
with dense computational grids in order to resolve the typical 
complex flow structure present in the jet entry into a cross- 
flow media. However due to the modeling assumptions and 
experimental uncertainty, variations in the numerically 
calculated results as compared to the measured data have 
been observed. It is therefore necessary to quantify the 
difference between the measured and numerically calculated 
mixedness results for a jet in a crossflow before further 


utilizing the analytical model as a configuration screening 
tool. In this investigation, the mixing effectiveness of a jet 
in a crossflow is numerically analyzed using a 3-D analytical 
code. Eighteen different configurations were analyzed 
which included a) circular orifices ' at three jet-to-manstream 
momentum-flux ratio (J) conditions, and b) slots with aspect 
ratios L/W of 4 and 8 with different slant angles between 0° 
and 90°. and at varying J conditions. The procedure used to 
model the geometric configurations was similar to that used 
by Oechsle. Mongia. and Holdeman. 1992. In this study 
however, the numerically calculated temperature field and 
temperature non-uniformity results were compared with 
measured data obtained by Hatch et al. (1992), and the 
difference has been quantified by measuring the variation of 
the temperature distribution with respect to the equilibrium 
temperature. 

Mathematical Model 

A production 3-D combustor code (COM-3D). solves the 
turbulent reacting flow transport equations using the SIMPLE 
algorithm of Patankar and Spalding (Patankar. 1980). This 
program simulates turbulence by the two-equation k-e model 
(Launder and Spalding, 1974), and combustion following 
vaporization is determined by a four-step chemical reaction 
model based on Arrhenius and eddy breakup concepts (Mon- 
gia and Reynolds. 1979). The transport equations for all 
dependent variables are of the following form as shown in 
Eq-1: 

div[ pu £, - (Meff/Pr) grad(^)] = S4 ( 1 ) 

where p r is the mixture density, u is the velocity. is the 
effective turbulent viscosity, P r is the effective 
Prandtl/Schmidt number, and S^ is the source term for the 
variable £. The following variables are computed by the 
three-dimensional code: 1) axial, radial, and swirl velocity 
components: 2) specific enthalpy and temperature; 3) turbu- 
lence kinetic energy and dissipation rate; 4) unburned fuel 
and composite fuel fraction; and 5) fuel spray trajectory and 
evaporation rate. 

The computational effort is significantly reduced by model- 
ing a sector of the mixing section comprising a single orifice. 
Therefore, the shape of the sector was dependent on the 
number of orifices equally spaced in the circumferential 
direction. It should be noted that there are 8 orifices in all 
configurations, so the computational domain is a 45 degree 
sector. Periodic boundary conditions were applied in the 
circumferential direction. These conditions were necessary 
for the mixing configurations generating swirl during the 
mixing process generally characteristic of the slanted slot 
orifice configurations. No-Slip and adiabatic boundary 
conditions were applied at the outer wail defining the inside 
wall of the mixing section. The computation at the center 
axis was extrapolated based on the nearest surrounding 
values due to its singularity. No-gradient and slip boundary 
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conditions were applied at the center axis. Fully developed 
profiles (no axial gradients) were assumed at the exit bound- 
ary condition. 


Geometric Configuration 

In this study, the mixing section was modeled as a constant 
diameter cylindrical "can" with a single row of equally 
spaced orifices. The outer wall diameter is 76 mm (3 inch) 
and the axial length of the mixing section extended from 
x/R=-1.5 to x/R=6 where x is referenced from the leading 
edge of the orifice. Sufficient axial distance was provided 
both upstream and downstream of the orifice to prevent any 
entry or exit effects from modifying the flow structure in the 
computational domain of interest which is 0<x/R<2. The 
computational grid of the domain typically comprised of 
50,000 to 80.000 control volumes. The grid was typically 
denser near the orifice and near the outer wall to resolve the 
high velocity and temperature gradients resulting from the 
inlet of the crossflow jet. An orthogonal view of the compu- 
tational grid is shown in Figure- 1. The grid was also config- 
ured to allow smooth progressive volume change between 
adjacent control volumes to help speed up the convergence 
of the solution. 

The geometric configurations of the jet orifices is also shown 
in Figure- 1. Two different general shapes of orifices were 
modeled: round holes and slanted slots of 4:1 and 8:1 aspect 
ratios, and a total matrix of 18 different configurations were 
analyzed. Also, it should be noted that the 90° sector shown 
in the following temperature plots always includes orifice 
centers at 22.5° and 67.5°. The description of the configura- 
tions is tabulated in Table- 1 and shown graphically in 
Figure-2. Twelve of the configurations shown in Table- 1 
duplicate the initial conditions reported in Hatch. Sowa. 
Samuelsen. and Holdeman, 1992). The configurations for 
which the measured data are not reported by Hatch et al. are 
configurations 4. 5. 11. 12. 14. and 15. The blockage param- 
eter tabulated in Table 1 is defined as the circumferential 
projection of the orifice divided by the spacing between the 
orifice centers. Similarly the orifice trailing edge (T.E.) x/R 
is defined as the axial projection of the orifice divided by the 
mixer radius R. As shown in the aforementioned experimen- 
tal study, the matrix of analyzed configurations is such that 
for a specific orifice shape, the increase in J does not change 
the MR. This was accomplished by reducing the orifice 
effective area (ACd) to increase the J value while keeping the 
MR constant (see Table- 1). The model of the orifice perim- 
eter was performed by a stair-stepping approximation since 
COM-3D is not a body conforming code. Therefore, to 
increase the modeling accuracy, the orifice was typically 
defined with between 100 and 200 control volumes. The 
orifice ACd was modeled to within less than 1 % deviation. 

The mixing section was modeled at atmospheric pressure. 
The mainstream flow and jet flow conditions are show in 
Table- 1. The MR. mainstream flowrate, and jet temperature 


were maintained constant throughout the parametric study. 
The variation in J was chosen as a representative range for 
typical gas turbine combustion systems. The mainstream 
flow at the inlet of the modeled mixing region was character- 
ized by a uniform temperature and axial plug flow velocity 
profiles in the radial-circumferential plane. The air jet How 
was characterized by a radial, uniform flow across the orifice 
effective area. The assumption of uniform mass injcc- 
tion/area is applied in the mathematical model in all the 
analyzed configurations. The turbulence kinetic energy of 
the mainstream and jet flows were 0.3% of the square of the 
mean velocities. The turbulence length scales of the main- 
stream flow were 2% of the can diameter, and the turbulence 
length scale of the jet was of the order of the orifice diameter. 
The inlet conditions for all orifices in the mixing section 
were equal to create a symmetrical input condition about the 
circumferential direction which was necessary for the sector 
analysis. The operating pressure and temperature range was 
chosen to duplicate the experimental setup. 

A typical numerical solution took about 250 iterations for full 
convergence with overall mass flow residuals of 0.05% of the 
total mixing section mass flowrate. All solutions were 
obtained using the Cray Y-MP and a typically converged 
solution took about 1 to 1.5 hrs of CPU time. 

Results and Discussion 

The mixing performance for all configurations analyzed in 
this study were ultimately evaluated at x/R =1. It is however 
recognized that the mixing performance throughout the entire 
mixing section 0<x/R<l should also be considered since 
complex structures are present in the flow field especially 
near the entry of the jet. The numerically calculated results 
were obtained using a grid with about 1200 computational 
nodes in each radial-tangential plane (45° sector). Similarly 
the grid size of the experimental data is 50 points per 
radial-tangential plane (90° sector) therefore in order to 
appropriately weigh the numerically calculated results, only 
the calculated results at the measured locations were used for 
the comparison. 

Two different methods were used to present the numerical 
results and compare them with their respective measured 
data: 

1) the results were analyzed qualitatively by visual obser- 
vation of the temperature and velocity field solutions. 
The normalized temperature plots presented in this paper 
describe a domain from x/R= 0.08 to x/R=l. In all the 
temperature plots, the air flows from left to right. The 
temperature results are presented as normalized values 
with respect to the overall differential between the main- 
stream flow and the jet flow inlet temperatures. This 
normalized temperature is defined in Eq-2. 
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f = JVl&l (2 ) 

T - T 
1 main A jet 

The value of f varies from 0 to 1, where 0 is the value of 
the unmixed jet and 1 is the value of the mainstream 
flow. Note that f = 1-0, where 0 is as defined previously 
(Holdeman, 1991 ) and used elsewhere also. The value of 
the equilibrium temperature was calculated in a purely 
adiabatic system at any location downstream of the jet 
injection. The equilibrium temperature was also normal- 
ized in the same manner as the local temperature. Five 
different planes are displayed in each plot corresponding 
to the planes in the radial-tangential direction at x/R of 
0.08, 0.25, 0.5, 0.75, and 1. Two adjacent sectors are 
plotted in each plane to show the continuity in the region 
between the orifices. Note that these planes were specifi- 
cally chosen since these are the actual measured planes. 
Each of the normalized temperature contour plots of the 
configurations with corresponding measurements include 
the following: 

a) the complete numerically calculated normalized 

temperature distribution 

b) the numerically calculated results at the measured 

locations 

c) plot of the measured results (Hatch, Sowa, Sam- 

uelsen, and Holdeman 1992.) 

2) The mixing performance of the mixing section was 
also evaluated by using statistical area weighted devia- 
tion. The smallest deviation indicates the best mixing 
configuration. Three different parameters (mixing devia- 
tions) AMIX. AHOT. and ACOLD are described in Eq-3 
through Eq-5. Note that AMIX also equals square root of 
the sum of the squared values of both AHOT and 
ACOLD. It is also important to note that the definitions 
in Eq-3 through Eq-5 do not correct the mixing 
non-uniformity for the bias introduced in the region of 
the section where the air is actively being injected 
through the orifice. This is accepted since the compari- 
son between the measured and numerically calculated 
results is only performed at x/R=l which is downstream 
of the trailing edge of the orifice. The numerically calcu- 
lated and measured area-weighted non-uniformity results 
are shown in Table-2. 


AMIX 


AHOT 



1/2 

(3) 


1/2 

for Tjk > T*, (4) 


ACOLD = 



Tjk - T«q 


1/2 

for T,k < T„ (5) 


In addition, the area weighted statistical deviations were also 
compared with a mass flow weighted parameters as defined 
in Eq-6 through Eq-8. The MMIX value also equals the 
square root of the sum of the squared values of both MHOT 
and MCOLD. Note that these parameters are essentially the 
same as the area weighted parameters with the density and 
velocity weighting terms added. 


MMIX = 


MHOT = 


MCOLD = 



Tjk - T«q 


Tjk - 

Tnui„ 




2 


2 


1/2 


1/2 

for T^ >T„ 


for Tj* <T eq 


( 6 ) 

(7) 

( 8 ) 


The results in this paper are presented in the following man- 
ner: 

a) effect of momentum flux ratio (J) on mixing. 

b) effect of orifice configuration on mixing, 

c) effect of the slot slant angle (<|)) on mixing, and 

d) comparison of the mixing trends obtained with the 
numerically calculated and measured results. 


Effect of Momentum Flux Ratio (J1 

Round holes. The mixing uniformity of all three round 
hole configurations are shown in Figures-3, 4, 5. and 6a for 
increasing J of 25. 52. and 80. It is important to note that the 
calculated values at the measured locations are essential to 
make a valid comparison between the numerically calculated 
and measured results. Since each quadrant in the measured 
data contains 50 points per plane, it will not resolve the large 
gradients in the temperature distribution usually present near 
the jet entrance. Therefore, in order to compare the measured 
and the numerically calculated results, the calculated data 
was linearly interpolated at the actual measured locations. 
The comparison of the experimentally calculated and mea- 
sured results for these configurations are also shown in 
Figure-7a. The lowest possible AMIX. AHOT, and ACOLD 
yields the best temeprature uniformity. The bar charts shown 
in Figure-7 depict the absolute value of the difference 
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between the mixing non-uniformity of the measured and 
numerically modeled statistical parameters. It is evident from 
the temperature plots that the increase of J beyond the value 
of 25 produces over-penetration. The mixing non-uniformity 
values AMIX, and AHOT both indicate a deteriorated perfor- 
mance with both the measured and calculated results which is 
in agreement with the expected results. However, the rela- 
tive comparison of the numerical and experimental results 
indicate that a better comparison is evident with increasing J 
as shown in Figure-7a (showing the decreasing relative 
difference between both values with increasing J). As previ- 
ously reported in Hatch et al. (1992), the increase in J pro- 
duced over-penetration which generally degraded the mixing 
uniformity. Similarly, the numerically calculated results also 
indicate some deterioration in mixing uniformity; however, 
this deterioration is smaller with increasing J (Figure-6a). 
There appears to be a slight improvement in the numerically 
calculated mixing uniformity results at J=80 due to the strong 
recirculation in the axial direction produced by the impinge- 
ment of opposing jets in the mixing section. The recircula- 
tion produces some additional jet-to-mainstream flow inter- 
action and additional mixing therefore reducing the effect hot 
mainstream air flow behind the wake of the jet. This effect 
appears to be diminished in the measured results since jet 
penetration was lower as compared to the calculated results 
at the same J; and therefore, the jet impingement recircula- 
tion is also smaller in the measured case. Note that the opti- 
mum jet penetration as calculated in Holdeman (1991) for a 
nominal J=25 is 9 holes/row as calculated in Eq-9 where n is 
the number of orifices and C=2.50. and these results seem to 
correlate well with the round jet results shown in Figure-3. 
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Slot aspect ratio L/W=4, with <b=45°. The mixing 
non-uniformity and normalized temperature profile plots for 
the 45° slanted slot with L/W=4 are shown in Figures-6b, 7b, 
8, 9, and 10. A significant difference between the measured 
and numerically modeled results are observed for the L/W=4 
slots as J is increased from 25 to 80 (Figure-7b). As J is 
increased, the swirl increases in the flow field. The swirl is a 
consequence of the jet flow acting like an effective turning 
vane which causes the mainstream flow to swirl in the cir- 
cumferential direction. Since there is no tangential velocity 
component to either the approach flow or jet. the induced 
swirl is actually a double counter-swirling flow such as to 
conserve the angular momentum of the entire system. The 
effect of this double recirculation shown in Oechsle, Mongia. 
and Holdeman (1992), can have a detrimental effect on the 
mixing uniformity in the numerical model since the main- 
stream flow becomes trapped in the wake of the jet and 
remains in the outer section of the mixer as the flow proceeds 
downstream. This creates a radially stratified temperature 
profile with hotter gas towards the outer wall of the mixing 
section. Due to the counter-swirling flow in the mixer down- 
stream of the orifice, the hot gas has little interaction with the 


cooler jet flow (only if the jet significantly penetrates into the 
inner recirculation core) since the stagnation region between 
the two counter-swirling flows produces an effective "bar- 
rier" to prevent significant mass flow transport in the radial 
direction. This effect appears to have a significant influence 
in the results shown in Figures-6b and 7b. As J is increased 
from 25 to 52. the radial temperature stratification becomes 
more significant in the numerically calculated results as 
compared to the measured results, and this appears to have a 
direct effect on the mixing non-uniformity since the hotter 
gas resides towards the outer section of the mixer. In the 
numerically calculated results, the increase in J from 25 to 80 
has a 10% decrease in the interface (r int ) between the 
counter-swirling flows as shown in Figure- 1 1 but the increase 
in jet penetration, as shown in Eq-9. is proportional to the Vj. 
and therefore jet penetration rather than r jm is more sensitive 
to J. It is therefore apparent that the numerically calculated 
results showing these temperature stratifications will deviate 
more compared to the measurements for orifice configura- 
tions with J=52 as shown in Figure-7b. As J is further 
increased to 80. the jet impingement recirculation produced 
by opposing jets aid the mixing the numerical results and 
therefore shows a slight improvement in mixing. 

Slot aspect ratio LAV=8, with 6=45°. The mixing unifor- 
mity and temperature profile plots for the 45° slanted slot 
with L/W=8 for increasing J from 25 to 80 are shown in 
Figures-6c. 7c, 12. 13. and 14. The penetration of the jet of 
the L/W=8 slot is relatively smaller as compared to the 
L/W=4 slots at the same J condition due to the long and 
narrow orifice configuration; and therefore, this is the only 
orifice configuration shape (of those analyzed herein) that 
approaches optimum penetration with increasing J. In addi- 
tion, the variation of jet penetration with J appears to be 
weaker as compared to the L/W=4 slots or circular holes. 
The numerically calculated results indicate a slight degrada- 
tion in the mixing non-uniformity with the increase in J 
(Figure-6c). Large differences in the comparison between 
the measured and calculated results are shown in Figure-7c. 
Of the cases analyzed numerically, none showed jet impinge- 
ment recirculation for the range of J of 25 to 80 and therefore 
no additional mixing benefit was derived as J is increased to 
80 as shown previously for the circular holes and L/W=4 
slots. Considering that the mixing development throughout 
the mixer downstream of x/R=0.5 is small, the deviation 
between the calculations and measurements will increase, 
especially if the measured results show significant mixing 
progress beyond x/R=0.5 as shown in Figures-12. 13, and 14. 


Effect of Orifice Configuration 


Considering that the design parameter J will most likely be 
constant for a given combustion liner design, it is therefore 
valuable to compare the mixing non-uniformity results for 
several orifice configurations at constant J. The round holes. 
L/W=4 slot, and L/W=8 slot orifice configurations are com- 
pared at J=25, 52, and 80 as shown in Figure- 15 for the 
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numerically calculated results. The temperature profiles 
corresponding to Figures- 15a and 16a are shown in 
Figures-3, 8. and 12. The temperature profiles corresponding 
to Figures-15b and 16b arc shown in Figurcs-4. 9. and 13. 
Similarly, the temperature profiles corresponding to 
Figures- 15c and 16c are shown in Figures-5. 10, and 14. 
These results almost exclusively show the L/W=8 slot as the 
best mixing configuration for all J values due to the shallow 
jet penetration and ability to mix well in the wake behind the 
orifice. It is important to note that the effect of jet 
under-penetration (produced by the L/W=8 slot) creates 
enhanced mixing towards the outer wall of the mixer at 
expense of a slightly hotter core (Figure- 12b) as compared to 
the round and L/W=4 slot results (Figures-3 and 8). The 
difference between the numerically calculated results and 
measured results show the large expected differences at 
<j>=45° as mentioned previously especially for J=52 as shown 
in Figure- 16. 


Effect of Slot Slant Angle (d>) 

The effect of the increase in slot slant angle on the mixing 
uniformity is shown in Figures- 17. 18, 19, 20. and 21 for the 
LAV =4 slot at slant angles of 0°, 22.5°, 45°. and 67.5° respec- 
tively for J=52. The difference in the mixing non-uniformity 
between the numerically calculated and measured results is 
shown in Figure-21c. As shown previously in Oechsle. 
Mongia, and Holdeman (1992). the increase in slot slant 
angle increases the induced flow counter-recirculation in the 
mixer which generally radially stratifies the temperature 
distribution (hotter towards the outer wall). This phenomena 
was mostly observed for the slant angle producing the stron- 
gest swirl (<t>=45°). It is also not surprising that the numeri- 
cally calculated mixing non-uniformity at <J)=45° is the 
highest as shown in Figure-21a. As <j> is increased to 67.5°, 
the jet penetration decreases sufficiently to allow the cold jet 
to mix with the hot gas trapped behind the orifice. The 
decrease in jet penetration with increased <J> beyond 45° is 
mainly caused by the increase in jet blockage, as mentioned 
previously in Oechsle. Mongia. and Holdeman (1992). The 
enhanced mixing for <p=67.5° is shown in Figure-21a. This 
effect appears to be a key element in producing a well mixed 
numerically calculated flow field. The measured flow field 
however, is generally well mixed for the slanted slot configu- 
rations with <J> from 22.5° to 67.5° and no significant change 
in performance was shown at J=52 (Figure-2 lb). 

The effect of <|> on the comparison of the numerically calcu- 
lated and measured mixing non-uniformity for J=25 and 80 
are not reported herein due to the lack of measured data 
reported, however, the numerical results were obtained at 
several (j) for the L/W=4 slot at J=25 and 80 (Figure-2). At 
J=25, the effect of mixing non-uniformity for a variation of <}> 
from 0 to 45° is shown in Figures-22 and 23a. It should be 
noted that at the low J condition, the 67.5° and 90° slots were 


not analyzed due to phisical constraints since the orifices will 
overlap each other. The results indicate small changes 
similar to those reported by Oechsle. Mongia. and Holdeman 
(1992). At J=80. the range in <|> from 0° to 90° (transversely 
oriented) for L/W=4 slots are shown in Figures-23c and 24 
respectively. The results are also similar to those reported in 
Oechsle. Mongia. and Holdeman ( 1992). The improvement 
in mixing for <|> above 45° is caused the same effect as 
explained for the J=52 results (Figure-23b). 


Comparison between the mixing trends for the numerically 
calculated and measured results. 

In general, the numerically calculated results follow the 
measured trends throughout the mixing section domain as 
seen in the Figures-25. 26. 27, and 28. It is evident from all 
the previously shown temperature plots that significant 
temperature gradients are still present in the numerically 
calculated results at x/R=l. This observation is expected, 
given the nature of the production elliptic code which has 
been observed and reported significantly in the past. 

Reiterating that lowest possible AMIX. AHOT. ACOLD. 
MMIX, MHOT, and MCOLD values depict best mixing, the 
results indicate similar trends for the measured as well as for 
the numerically calculated results in the development of the 
mixing flow field. However the temperature gradients in the 
calculations appear not to significantly change downstream 
of x/R=0.5 indicating very little mass or thermal transport in 
the radial and circumferential directions downstream of 
x/R=0.5. This was not the case with the measured results. 
The mixing non-uniformity parameters obtained with the 
numerically calculated data tire generally higher compared to 
the measured values which is due to the difference in the 
higher temperature gradients. This also indicates that the 
numerically modeled results are less mixed compared to the 
measured values. 

In addition to the area weighted results presented in 
Figures-25 through 28. the mass flow weighted parameters 
are also shown. The results indicate that the general trends 
are also comparable between the area and mass flow 
weighted parameters: however, significant deviations occur 
at locations x/R<0.8 which is in the domain of jet injection. 
These results are expected to be different since the area 
weighted results do not account for the increment in mass 
injection through the jet. Since the mass flow weighted 
results contain the axial velocity term, they reflect this 
difference. Note that the mass flow weighted deviations for 
the measured data are not reported herein since the velocity 
field was not measured by Hatch et al (1992). 

The best mixing configuration based on the numerically 
calculated results shown herein at x/R=l is the 67.5° slanted 
slot with LAV=4 at J=52 (Figures-20, and 28j. 28k, and 281). 
Similarly the best measured configuration with regard to 
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mixing uniformity at x/R=l is the slanted slot configuration 
with a L/W=4 slant angle range of <t>=45° at J=52 (Figures-9, 
28g, 28h, and 28i). Both numerically modeled and measured 
methods therefore show a good comparison as to the mixing 
uniformity at x/R=l. Note that the measured and numerically 
calculated overall performances were determined based on a 
numerical average of the AMIX, AHOT, and ACOLD param- 
eters for each specific configuration. 


Conclusions 

1) The area weighted parameters compare favorably with the 
mass flow weighted parameters in almost all numerically 
analyzed configurations at x/R=l. It is however noted that 
due to the definitions of these statistical parameters, differ- 
ences will be encountered in the region where the jet flow is 
injected (typically 0<x/R<0.8 for the cases analyzed herein). 

2) The difference between the numerically calculated and 
measured results increases with increasing circumferential 
counter-swirl (generally produced by slot slant angles near 
45° and high J) since the mixing development and tempera- 
ture uniformity improves in the measured results beyond 
x/R=0.5. The change in the numerically calculated tempera- 
ture gradients and mixing uniformity with increasing x/R 
appears to be small beyond x/R=0.5. 

3) The radial temperature stratification produced by the 
slanted slot configurations in the numerically calculated 
results is significant and therefore large variations in the 
mixing non-uniformity between these results and the mea- 
surements were encountered. In the numerically calculated 
results, optimum penetration was therefore obtained with 
relatively shallow jets which became entrained behind the 
orifice and mix with the entrained hot flow. 

4) The numerically calculated results indicate a much larger 
variation in mixing non-uniformity with a change in slot slant 
angle as compared to the measured results. For J=25, the jet 
penetration is shallow enough to mix cold jet flow with hot 
mainstream gas behind the jet near the outer wall therefore 
producing increased mixing and lower temperature stratifica- 
tion. As J is further increased, the jet penetrates beyond the 
recirculation stagnation therefore trapping the colder jet flow 
towards the inner section of the mixer and producing high 
temperature stratification towards the outer wall of the mixer. 

5) The increase in J. creating an over-penetrating jet situa- 
tion, slightly improves the mixing uniformity due to addi- 
tional mixing produced upstream of the orifice due to the jet 
impingement recirculation which allows cold air to be 
entrained behind the orifice and mix with the hot mainstream 
flow. 


6) The general trends shown in the numerically calculated 
results are similar to the measured results and the best overall 
mixing configurations are comparable; however, the com- 
parison between the numerically calculated and measured 
results must be made at the same locations to equally bias 
both results for an accurate comparison. 
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Orifice aspect ratio (LAV) 
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Orifice width 
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Orifice T.E. x/R 
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OVERALL 


J 


MR (jet/main) 


DR Get/main) 


VR ( iet/main 


AR (jet/main) 
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MAINSTREAM 


P main 


T main 


U main 


m main 

Mixer Diameter 


JET 


T jet 
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m iet 


ACd iet/row 


I Orifice Cd 

Number of orifices 


Orifice aspect ratio (L/W) 

Slant angle 

Orifice length 


Orifice width 

Blockage 

Orifice T.E. x/R 


99.3 
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28.1 

2.20 
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Table- 1. Overall description of the operating conditions and geometric dimensions 
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mixing non-uniformity results for all the analyzed configurations at x/R=l 




























































































































































Figure- 1 Depiction of (a) geometric configuration of the mixing 
section; (b) computational grid 
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Figure-2. Matrix of analyzed configurations 



b. Calculated results at the measured locations 
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Figure-3. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 1, J=26.7 MR- 2.2, DR=1.26, and T e =0.313 

Round jet , 8 orifices/row 
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a. Numerically calculated results 



b. Calculated results at the measured locations 
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Figure-4. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 2, J=55.4 , MR-2.2, DR= 1.26, and T t =0.313 

Round jet, 8 orifices/row 
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a. Numerically calculated results 



x/R 


b. Calculated results at the measured locations 
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Figure-5. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 3, J-84.2 , MR- 2.2, DR- 1.26, and T. =0.313 

Round jet, 8 orifices/row 
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MIXING NON-UNIFORMITY 



Nominal Momentum flux ratio (J) 
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Figure-6. Effect of the increase in J on the numerically calculated mixing 
non-uniformity of round holes and slanted slots with aspect ratios of 4 and 8. 
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Figure-7. Effect of the increase in J and slant angle on the difference between the 

numerically calculated and measured results for different orifice configurations, at x/R=l. 


22 


1.00 


a. Numerically calculated results 
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b. Calculated results at the measured locations 
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Figure-8. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 6 , J=30.5, MR=2.2, DR= 1.26 , and T =0.313 

45 0 slot , U\V=4 , 8 orifices/row 
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a. Numerically calculated results 
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b. Calculated results at the measured locations 
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Figure-9. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 9, J=57.7 , MR-2.2, DR= 1.26 , and T e =0.313 

45° slot, L/\\=4, 8 orifices/row 

25 na . iHiiHuomm m**# 




BLANK KOT FILMED 




a. Numerically calculated results 
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Figure-10. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 13, J-93.0, MR=2.2, DR=1.26, and T e =0.313 

45° slot, UW-4, 8 orifices/row 
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a. Velocity distribution in the 
radial-tangential plane at x/R=l 
for configuration #18, 45 deg 
slots, L/W=8 
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Figure- 1 1 Effect of J on counter-swirl interface radius 
for the numerically calculated resutls 
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a. Numerically calculated results 



b. Calculated results at the measured locations 





Figure-12. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 16, J-26. 7 , MR=2.2, DR-1.26, and T t( = 0.313 

45° slot, L/W= 8, 8 orifices/row 
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a. Numerically calculated results 



Figure-13. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration #17, J=50.9, MR- 2.2, DR=1.26, and T e = 0.313 

45° slot, L/W= 8, 8 orifices/row 
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a. Numerically calculated results 



Figure-14. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 18, J=88.5, MR=2.2, DR=1.26, and T eq =0.313 

45° slot, L/W=8, 8 orifices/row 
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Figure- 15. Effect of orifice configuration and slot aspect ratio (LAV) on the 
numerically calculated mixing non-uniformity for varying J. 
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Figure- 16. Effect of orifice configuration and slant slant angle on the difference 
between the numerically calculated and measured results at x/R=l. 
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a. Numerically calculated results 
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b. Calculated results at the measured locations 




Figure-17. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration # 7, J=5l.l , MR-2.2, DR= 1.26 , and T c =0.313 

Aligned slot , L/W=4 , 8 orifices/row 
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a. Numerically calculated results 
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b. Calculated results at the measured locations 
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Figure-18. Comparison of the normalized temperature for the numerically calculated 


and measured results for configuration #8, J=53.2 , MR=2.2 , DR= 1.26, and T e =0.313 


22.5° slot, U\\=4 , 8 orifices/row 
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a. Numerically calculated results 



x/R 


b. Calculated results at the measured locations 
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Figure-19. Comparison of the normalized temperature for the numerically calculated 
atid measured results for configuration # 9 , J=57.7, MR=2.2, DR=1.26, and T e = 0.313 

45° slot, UW= 4, 8 orifices/row * 
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a. Numerically calculated results 



Figure-20. Comparison of the normalized temperature for the numerically calculated 
and measured results for configuration #10, J-59.9, MR=2.2, DR- 1.26, and T eq - 0.313 

67.5 ° slot, L/W=4, 8 orifices/row 
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Figure-21. Effect of slot slant angle (<j> ) on the difference between the 
calculated and measured results for the L/W=4 slots at J=52 
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b. 22.5° slot , L/W=4, 8 orifices/row, J-28.0 
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Figure-22. Normalized temperature distribution for the numerically calculated 
results for MR=2.2, DR= 1.26, and T = 0.313 
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MIXING NON-UNIFORMITY 



Figure-23. Effect of the slot slant angle (<|> ) on the numerically calculated 
mixing results for the L/W=4 slots. 
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Figure-24. Normalized temperature distribution for the numerically calculated 
results for MR-2.2, DR= 1.26 , and T ,=0.31 3 
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Figure-24 (continued) 
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Conf # 1 , round holes, J=26.7, MR=2.2, DR=1.26 (a, b, c) 
Conf # 2 , round holes, J=55.4. MR=2.2. DR=1.26 (d, e, f) 
Conf # 3 , round holes, J=X4.2, MR=2.2, DR=1.26 (g, h, i) 
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Figure-25. Comparison of the numerically calculated and measured results 

for round holes. 
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Conf # 6, 45° slot, L/W=4, J=30.5, MR=2.2, DR=1.26, (a,b,c) 
Conf # 9, 45° slot, L/W=4, J=57.7, MR=2.2, DR=1.26, (d,e,f) 
Conf# 13, 45° slot, L/W=4, J=93.0, MR=2.2, DR=1.26, (g,h,i) 
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Figure-26. Comparison of the numerically calculated and measured results 
for the 45° slanted slot, L/W=4 . 
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Mixing non-uniformity 










Conf # 16, 45° slot, L/W=8, J=28.1, MR=2.2, DR=1.26, (a,b,c) 
Conf # 17, 45° slot, L/W=8, J=50.9, MR=2.2, DR=1.26, (d,e,f) 
Conf# 18, 45° slot, L/W=8, J=88.5, MR=2.2, DR=1.26, (g,h,i) 
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Figure-27. Comparison of the numerically calculated and measured results 
for the 45° slanted slot, L/W=8 . 
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Conf # 7, 0° (aligned) slot, L/W=4, J=51.1, MR=2.2, DR=1.26, (a,b,c) 
Conf # 8, 22.5° slot, L/W=4, J=53.2, MR=2.2, DR=1.26, (d,e,f) 

Conf # 9, 45° slot, L/W=4, J=57.7, MR=2.2, DR=1.26, (g,h,i) 

Conf# 10, 67.5° slot, L/W=4, J=59.9, MR=2.2, DR=1.26, (g,h,i) 
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Figure-28. Comparison of the numerically calculated and measured results 
L/W=4 slots at a nominal momentum flux ratio J=52. 
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